Collective effects in intra-cellular molecular motor transport: 
coordination, cooperation and competetion * 
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Molecular motors do not work in isolation in-vivo. We highlight some of the coordinations, coop- 
erations and competitions that determine the collective properties of molecular motors in eukaryotic 
cells. In the context of traffic-like movement of motors on a track, we emphasize the importance of 
single-motor bio-chemical cycle and enzymatic activity on their collective spatio-temporal organi- 
sation. Our modelling strategy is based on a synthesis- the same model describes the single-motor 
mechano-chemistry at sufficiently low densities whereas at higher densities it accounts for the col- 
lective flow properties and the density profiles of the motors. We consider two specific examples, 
namely, traffic of single-headed kinesin motors KIF1A on a microtubule track and ribosome traffic 
on a messenger RNA track. 
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I. INTRODUCTION 

Molecular motors 0, perform crucial functions at 
almost every stage in the life of a cell. The motor- 
dependent activities begin with DNA replication Q 
which is essential before cell division. A motor called 
DNA helicase 0, |B| unzips the double-stranded DNA 
so that another motor, called DNA polymerase Q can 
synthesize two copies of the DNA using the two single 
strands as templates. Protein synthesis one of the 
crucial activities within a cell, is accomplished by some 
other motors at different stages of operation, the most 
notable among these are RNA polymerase Q and ribo- 
somes 0. Normally, the proteins are synthesized near 
the cell center and need to be transported to appropri- 
ate locations some of which are far from the cell center. 
Besides, some cargoes are also transported from the cell 
periphery to the cell center. Most of these transporta- 
tions of cargoes (often in appropriately packaged form) 
in the cytoplasm are carried out by yet another set of 
molecular motors which are collectively referred to as cy- 
toskeletal motors [(|. 

The molecular motors at the cellular and subcellular 
levels can be grouped together is several different ways. 
One possible way of grouping these is based on the na- 
ture of the tracks on which the motors move- one group 
moves on filamentary proteins whereas the other moves 
on nucleic acid strands. Among the cytoskeletal motors, 
which walk on filamentary proteins, kinesin and dynein 
move on microtubules (MT) whereas myosins move on 
actin filaments. Some motors are processive, i.e., capable 
of long walk without getting detached from the track. Ki- 
nesins and dyneins are not only processive but can carry 



molecular cargo; therefore, these are functionally similar 
to "porters". 

Molecular motors do not work in isolation in-vivo. 
Wide range of biophysical phenomena observed at sub- 
cellular and cellular levels are manifestations of different 
types of collective processes, involving molecular motors, 
at several different levels of biological organisation ^(J . 
The nature of the collective effect depends on the situ- 
ation; it can be the coordination of different parts of a 
single motor or, on a larger scale, it can be a cooperation 
or competition between two different motors. These col- 
lective effects can manifest as spontaneous oscillations, 
dynamic instabilities, hysteresis, dynamical phase tran- 
sitions 0,0] and motor traffic jam In this paper, 
we first give examples of different types of coordination, 
cooperation and competition in molecular motor trans- 
port. Then, we focus on our works on molecular motor 
traffic on two different types of tracks, namely, on MT 
and on messenger RNA (mRNA). 



II. COORDINATION THROUGH ELASTIC 
COUPLING 

In this section we consider coordination of elastically 
coupled motor domains of a given motor as well as that of 
motors which are elastically coupled to a common back- 
bone. The structures and mechanisms of single cytoskele- 
tal motors have been discussed in this proceedings by Ray 
H3, and by Mallik and Gross [3 while those of helicase 
motors have been reviewed by Tuteja and Tuteja [Tfil ]. 



Coordination of different heads of oligomeric 
motors 



* Section IV Bl is based on original work l59t l67l carried out in col- 
laboration with K. Nishinari, Y. Okada, A. Schadschneider, P. 
Greulich and A. Garai. Section IVII is based on original work l7Sl 
done in collaboration with A. Basu. 



Most of the known common molecular motor proteins 
are dimeric or, in general, oligomeric 17]. It is generally 
believed that coordination of the ATPase activities of 
different motor domains is essential for the processivity 
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of a given motor. The mechanism of this coordination has 
been investiagted extensively over the last decade using 
many different techniques. 

It is now quite well established that the conventional 
double-headed kincsins follow a hand- over- hand mecha- 
nism (exactly similar to the steppings used by humans for 

walking) [lllHi^Eililii. Similar h and-over-hand 
mechanism is also believed to govern the movement of 
myosin-V, an unconventional processive myosin p3 . |25| . 
However, the corresponding mechanism for dimeric cyto- 
plasmic dynein is much more complex because the archi- 
tecture of its motor domain is very different from those 
of kinesin and myosin-V j2|| [2^ • 

Majority of the molecular motors that move on nucleic 
acid tracks (DNA or RNA) are also oligomeric. For exam- 
ple, most DNA helicases are either dimeric or hexameric 
|4(. Both "inchworm" and "hand-over-hand" (also called 
"rolling" ) mechanisms for the coordination of the motor 
domains of dimeric DNA helicases have been considered 
|28j. The corresponding mechanism for hexameric DNA 
helicases remains unclear; at least three different plausi- 
ble scenarios have been suggested. The ATPase activity 
of the motor domains in hexameric DNA helicases can 
run (i) in parallel, or (ii) in ordered sequential, or (hi) in 
random-sequential manner 29] . 



B. Collective dynamics of cytoskeletal motors 
bound to a common elastic backbone 

Consider a group of identical motors bound to an elas- 
tic backbone as shown in figQ] It has been shown [3(j 
that, even if each individual motor is non-processive, such 
a system of elastically coupled motors can move collec- 
tively on a filamentary track in a processive manner in 
one direction for a period of time and, then, sponta- 
neously reverse its direction of motion. Such spontaneous 
oscillations can account for the dynamics of axonemes, 
which are core constituents of eucaryotic cilia, as well 
as oscillatory motions of flight muscles of many insects 

MM- 



III. TWO SUPERFAMILIES OF PORTERS: 
LOAD-SHARING VERSUS TUG-OF-WAR 

Kincsins and dyneins are both processive (i.e., 
"porters"), but move in opposite directions along MT 
tracks. Often several members of the same superfam- 
ily together carry a single cargo. On the other hand, 
if members of both the superfamilies are adsorbed si- 
multaneously on the same cargo, they compete against 
each other trying to pull the cargo in opposite directions. 
Recent progress in theoretical modelling of these phe- 
nomena have been summarized in this proceedings by 
Lipowsky et al.|32T] 




FIG. 1: Schematic descri ptio n of the model of elastically cou- 
pled cytoskeletal motors |3C| : the horizontal curve represents 
the elastic backbone. 

A. Unidirectional transport: load-sharing by 
members of one superfamily 

Two types of in-vitro motility assays have been used 
for studying cytoskeletal molecular motors. In particu- 
lar, in the bead assay the filamentary tracks are fixed to 
a substrate and motors are attached to a micron-sized 
bead (usually made of glass or plastic). The movement 
of the bead in the presence of ATP is monitored using 
appropriate optical micropscopic methods. In such sit- 
uations, each bead is likely to be covered by N motors 
where, in general, N > 1. More than one motor is also 
used for transportation of large organelles in-vivo. This 
phenomenon of unidirectional cooperative cargo trans- 
port by more than one motor of the same superfamily 
has been studied theoretically [13 • The average walking 
distance increase with increasing N. In the face of oppo- 
sition by external load force, the force is shared by the 
N motors. 

Cooperative pulling of lipid bilayer membranes by ki- 
nesin motors can generate membrane tubes |34i l35l l36| . 
Such minimal systems are adequate to gain insight into 
the role of motor-mediated interaction between the cy- 
toskeleton and organelles in the formation of endoplasmic 
reticulum and Golgi membrane networks. 

In case of transport by cytoplasmic dynein motors, if 
more than one motor work together then collectively they 
can improve the performance as porter although individ- 
ually, while working in isolation, they often pause and 
experience slippage in presence of load force [3jj • 

B. Bidirectional transport: tug-of-war between two 
superfamiles 

In the preceeding section we have mentioned how the a 
single assembly of elastically coupled motors can sponta- 
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FIG. 2: Schematic description of the three possible mecha- 
nisms of bidirectional traffic. The three figures from top to 
bottom correspond to the three mechanisms (i), (ii) and (iii) 
explained in the text. The open circle and the rectangular box 
represent, respectively, a micron-size bead and a regulator of 
bidirectional transport. 



neously reverse its direction of motion on the same track. 
It is also well known that some motors reverse the di- 
rection of motion by switching over from one track to 
another which are oriented in anti-parallel fashion. In 
contrast to these types of reversal of direction of mo- 
tion, we consider in this section those reversals where the 
cargo uses a "tug-of-war" between kinesins and dyneins 
to execute bidirectional motion on the same MT track 
[3sL l3^ | . Several possible functional adva ntag es of bidi- 
rectional transport have been conjectured |38l l39| . 

Wide varieties of bidirectional cargoes have already 
been identified so far; these include organelles (for exam- 
ple, mitochondria) as well as secretaory vesicles and even 
viruses. If motion in one direction dominates overwhelm- 
ingly over the other, it becomes extremely difficult to 
identify the movement unambiguously as "bidirectional" 
because of the limitations of the spatial and temporal 
resolutions of the existing techniques of imaging. 

The main challenge in this context is to understand 
the mechanisms of this bidirectional transport and those 
which control the duration of unidirectional movement in 
between two successive reversals. This insight will also 
be utilized for therapeutic strategies. For example, the 
motor or the motor-cargo link may be targeted block- 
ing the virus that hijacks the motor transport system to 
travel towards the nucleus. On the other hand, a virus 
executing bidirectional movements can be turned away 
from the outskirts of the nucleus by tilting the balance 
in favour of the kinesins. 

At least three possible mechanisms of bidirectional 
transport have been postulated (see figEl- (i) One pos- 
sibility is that either only + end directed motors or only 
- end directed motors are attached to the cargo at any 
given instant of time. Reversal of the direction of move- 



ment of the cargo is observed when the attached mo- 
tors are replaced by motors of opposite polarity, (ii) The 
second possible mechanism is the closest to the real life 
"tug-of-war" ; the competion between the motors of op- 
posite polarity, which are simultaneously attached to the 
same cargo and tend to walk on the same filament gener- 
ates a net displacement in a direction that is decided by 
the stronger side, (iii) The third mechanism is based on 
the concept of regulation; although motors of opposite 
polarity are simultaneously attached to the cargo, only 
one type of motors are activated at a time for walking on 
the track. In this mechanism, the reversal of the cargo 
movement is caused by the regulator when it disengages 
one type of motor and engages motors of the opposite po- 
larity. For experimentalists, it is a challenge not only to 
identify the regulator, if such a regulator exists, but also 
to identify the mechanism used by the regulator to act 
as a switch for causing the reversal of cargo movement. 
Dynactin has been identified as a p ossible candidate for 
the role of such a regulator 0> l4l| . 



IV. TRANSFER FROM MT NETWORKS TO 
ACTIN NETWORKS: PARK-AND-RIDE 
TRANSPORT SYSTEM 

The networks of MTs and actin filaments are not dis- 
connected. The cytoskeleton is MT-rich near the cell cen- 
ter whereas dense actin filaments dominate the cytoskele- 
ton near the cell periphery. In recent years functional co- 
operation between the MT-based and actin-based trans- 
port systems have been discovered 0> EE EH ■ 

Several organisms are known to have an intrinsic abil- 
ity to manipulate their skin colour. For example, in 
fish and frog this change of skin colour is caused by 
the appropriate movement of pigment-containing vesicles 
called melanosome within a special type of cells called 
melanophore. When stimulated, melanosomes disperse 
throughout the cell causing darkening of the skin colour. 
A second stimulus can reverse the process by aggrega- 
tion of the melanosomes near the cell center. During dis- 
persion melanosomes are transported, some part of their 
way towards the cell periphery, by kinesin motors. Then, 
melanosomes switch from the MT-based kinesin motors 
to actin-based myosin motor (unconventional myosin- 
V which is a processive motor) for further distribution 
around the cell periphery. During the reverse process, 
the melanosomes begin their journey on the actin-based 
transportation network and, somewhere along the way, 
switch over to the MT-based network for continuing their 
onward journey towards the cell center. This switching of 
the transport system is similar to park-and-ride system- 
one drives a car from the suburban areas to the near- 
est station of the urban high-speed mass transit system 
(railways or metro trains) for reaching the central part 
of a city . 

Cooperation of MT-based and actin-based transporta- 
tion networks also play crucial roles in neuronal trans- 



port. Actin filaments help in the transport of the cargo 
to bridge the gaps between MT filaments |42j]. The link- 
ers of the two networks and the mechanisms of regulation 
and control of the proper switchover of cargo from one 
system to the other have received some attention in the 
last few years Efil l47j . 



V. COLLECTIVE TRANSPORT BY KINESINS: 
MOLECULAR MOTOR TRAFFIC JAM 

Most of the multi-motor phenomena we have consid- 
ered in the preceeding section are restricted to sufficiently 
low densities where direct interaction of the cargoes did 
not occur. As the cargoes are always much bigger than 
the motors (in- vitro as well as in-vivo), direct steric in- 
teractions of the cargoes become significant when several 
cargoes are carried by sufficiently dense population of 
motors along the same track. Such situations are remi- 
niscent of vehicular traffic where mutual hindrance of the 
vehicles cause traffic jam at sufficiently high densities. In 
analogy with vehicular traffic, we shall refer to the collec- 
tive movement of molecular motors along a filamentary 
track as "molecular motor traffic"; we shall explore the 
possibility of molecular motor traffic jam and its possible 
functional implications. 



A. TASEP-like minimal models of molecular motor 

traffic 

Most of the minimal theoretical models of interacting 
molecular motors 0> HH> EH EH E3 utilize the similari- 
ties between molecular motor traffic on MT and vehicular 
traffic on highways [5^ both of which can be modelled 
by appr opriate extensions of driven diffusive lattice gases 
[541 l55| . In such models the motor is represented by a 
"self-propelled" particle and its dynamics is formulated 
as an appropriate extension of the totally asymmetric 
simple exclusion process (TASEP) which is one of the 
simplest models of interacting driven particles. Many as- 
pects TASEP and its extentions as well as generalizations 
have been discussed in detail in the articles by Stinch- 
combe |56|. Kolomeisky |57l ] and Barma |58| in this pro- 
ceedings. The novel feature of these models, which distin- 
guish these from TASEP are the Langmuir-like kinetics 
of adsorption and desorption of the motors. 

Let us consider the model suggested by Parmeggiani et 
al.|4^. In this model (see figOJ) the self-propelled particle 
can hope from one lattice site to the next immediately in 
front of it, provided the target site is empty, with hopping 
probability q per unit time. Moreover, a particle gets at- 
tached to an empty site with probability A per unit time 
time, whereas a particle attached to a site can get de- 
tached with probability D per unit time. Furthermore, 
the attachment probability at the entrance and the de- 
tachment probability at the exit are denoted by a and (3, 
respectively, per unit time. In spite of the extreme sim- 
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FIG. 3: Schematic description of the model of cytoskeletal 
motor traffic developed by Parmeggiani et al. |49(l 



M = Motor 
T = Track 




M < ; ^ M.NTP 



NTP = 




FIG. 4: Biochemical cycle in the absence and presence of the 
track. The symbols above (or below) the arrows are the rate 
constants for to the chemical reactions represented by the 
corresponding arrows. 



plicity of this model, it predicted a novel jammed phase 
on the a— (3 phase diagram. The progress made so far has 
been reviewed in this proceedings by Lipowsky et al.[32j|. 



B. KIF1A traffic: effects of biochemical cycle 

In reality, a molecular motor is an enzyme that hydrol- 
yses ATP. In the presence of the track, usually the AT- 
Pase activity of the motor increases by one or two orders 
of magnitude; moreover, the mechanical movement of the 
motor is coupled to its biochemical cycle that includes its 
enzymatic activity (see figEJ. Models developed for de- 
scribing the mechanisms of operation of a single motor of 
a specific type incorporate not only the mechanical states 
of the motor but also its chemical states in each biochem- 
ical cycle. On the other hand, the models of interacting 
motors in traffic, which we mentioned in the preceeding 
subsection, do not explicitly take into account the dis- 
tinct chemical states of each motor during its biochemical 
cycle. Therefore, one common theme of our recent works 
on molecular motors has been to achieve a synthesis- we 
incorporate the essential steps of the biochemical cycle, 
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in addition to the mutual interactions of the motors, to 
develope models for their collective traffic-like dynamics. 
In the low-density limit, the model describes single- motor 
dynamics while the same model at higher densities pre- 
dicts the collective spatio-temporal organization of the 
motors. 

In a recent paper (5^ we considered specifically the 
single-headed kinesin motor, KIF1A |6Ct l6ll l62| ; the 
movement of a single KIF1A motor had already been 
modelled with a Brownian ratchet mechanism |63l l64l . 
In contrast to the earlier models 0, 14^. IHol l5lL l52| 
of molecular motor traffic, which take into account only 
the mutual interactions of the motors, our model (from 
now referred to as the NOSC model) explicitly incorpo- 
rates also the Brownian ratchet mechanism of individual 
KIF1A motors, including its biochemical cycle that in- 
volves ATP hydrolysis. 

The biochemical cycle of a single-headed kinesin mo- 
tor KIF1A can be described by the four states shown in 
Fig.[51 [60Ll62| : bare kinesin (K), kinesin bound with ATP 
(K.ATP), kinesin bound with the products of hydrolysis, 
i.e., adenosine diphosphate and phosphate (K.ADP.P), 
and, finally, kinesin bound with ADP (K.ADP) after re- 
leasing phosphate. 

Through a series of in-vitro experiments, Okada, Hi- 
rokawa and co-workers established that 

(i) each KIF1A motor, while weakly bound to a MT, re- 
mains tethered to the MT filament by the electrostatic 
attraction between the positively charged if -loop of the 
motor and the negatively charged i?-hook of the MT fil- 
ament. 

(ii) In the weakly bound state, a KIF1A cannot wan- 
der far away from the MT, but can execute (essentially 
one-dimensional) diffusive motion parallel to the MT fil- 
ament. However, in the strongly bound state, the KIF1A 
motor remains strongly bound to a motor binding site on 
the MT. 

(iii) A transition from the strongly bound state to the 
weakly bound state is caused by the hydrolysis of the 
bound ATP molecule. After releasing all the products of 
hydrolysis (i.e., ADP and phosphate), the motor again 
binds strongly with the nearest binding site on the MT. 

In the NOSC model a single protofilament of MT is 
modelled by a one-dimensional lattice of L sites each of 
which corresponds to one KIFlA-binding site on the MT; 
the lattice spacing is equivalent to 8 nm which is the sep- 
aration between the successive binding sites on a MT 
0. Each kinesin is represented by a particle with two 
possible internal states, labelled by the indices 1 and 2, 
in which it binds strongly and weakly, respectively, to 
the MT track. Attachment of a motor to the MT oc- 
curs stochastically whenever a binding site on the latter 
is empty. Such a two-state approximation to the full se- 
quence of biochemical states (conformations) has been 
successfully exploited also for conventional two-headed 
kinesin and unconventional my osin- V both of which are 
known to be processive [65, 66] . 

For the dynamical evolution of the system, one of the 



K.ADP.P 



K.ATP 




FIG. 5: The biochemical cycle of a single-headed kinesin mo- 
tor KIF1A. 



L sites is picked up randomly and updated according to 
the rules given below together with the corresponding 
probabilities (Fig.|f)J|: 



Attachment : 
Detachment 
Hydrolysis : 

Ratchet : 





1 • 

1 - 

2 
20 



Brownian motion 



-> 1 with uj a dt 

* with LOddt 

2 with ojhdt 

-> 1 with Lu s dt 
-> 01 with ujjdt 

20 -> 02 with Lu b dt 
02 -> 20 with w b dt 



(1) 
(2) 
(3) 

(4) 
(5) 



The physical processes captured by the rate constants 
Wf and w s can be understood as follows by analyzing 
the Brownian ratchet mechanism illustrated in figEl For 
the sake of simplicity, let us imagine that the potential 
seen by the motor periodically oscillates between the saw- 
tooth shape and the flat shape shown in figEl When the 
sawtooth form remains "on" for some time, the particle 
settles at the bottom of a well. Then, when the potential 
is switched "off" , the probability distribution of the posi- 
tion of the particle is given by a delta function which, be- 
cause of free diffusion in the absence of any force, begins 
to spread. After some time the Gaussian profile spreads 
to such an extent that it has some overlap also with the 
well in front, in addition to the overlap it has with the 
original well. At that stage, when the sawtooth potentil 
is again switched on, there is a non- vanishing probability 
that the particle will find itself in the well in front; this 
probability is proportional to the area of that part of the 
Gaussian profile which overlaps with the potential well in 
front and is accounted for in our model by the parameter 
W / . There is also significant probability that the particle 
will fall back into the original well; this is captured in 
our model by the parameter w s . 

Let us denote the probabilities of finding a KIF1A 
molecule in the states 1 and 2 at the lattice site i at time 
t by the symbols and hi, respectively. In mean-field 
approximation the master equations for the dynamics of 
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Assuming periodic boundary conditions, the solutions 
(ri,hi) — (r,h) of the mean- field equations J7J in the 
steady-state are found to be 



-n h - n s - (n s - i)k + Vd 

2K(1 + K) ' 

n h + si s + (n. + 1) k - Vd 

2K 

where K — ojd/u a , &h = ^h/^ft = ^s/^fj and 



(8) 
(9) 



FIG. 6: Schematic description of the Brownian-ratchet. 



i-1 



i+1 



d = 4n s K(i + k) + (fi ft + n s + (n s - i)K) . (io) 

The probability of finding an empty binding site on a 
MT is Kr as the stationary solution satisfies the equation 
r + h + Kr = 1. The corresponding flux is given by |67j 
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UJ a UJ h Uf 



Uf(u a + UJ d ) + U a (iJ B + UJ h ) + ^d^s 

u h (l + K) 

(i + K)[(i + k) + (n s + n h ) + n s K] ' 



(11) 



We have also computed the average density profile of 
the motors along a MT track with open boundary con- 
ditions. For a given u> a , the density of motors in state 2 
away from the edges of the MT exceeds that of those in 
state 1 as u>h increases beyond a certain value. 



FIG. 7: A schematic representation of the NOSC model of 
KIF1A traffic. The integer i denotes the discrete equispaced 
binding sites for kinesin motors on the MT track. The sym- 
bols 1 and 2 within circles at a given site i represent the 
"chemical" (conformational) states of a kinesin at the same 
spatial location. 



motors in the bulk of the system are given by |5£ 




= w a (l - n - hi) - oJhn - ojdn + oj s hi 

+cjfhi-i(l -n - hi), (6) 

^ = -cj s hi + u) h ri - Ufhi(l - r i+ i - h i+ i) 

-uj b hi(2 - n+i - h i+ i - n-x - hi-x) 
+uj b (hi-i +h i+ i)(l - n - hi). (7) 

The corresponding equations for the boundaries, which 
depend on the rate constants of entry and exit at the 
two ends of the MT, are similar and will be presented 
elsewhere [67|. 

The sing le KIFIA properties are reproduced by this 
model |59| in the appropriate low-density limit. For 
example, v, the mean speed of the kinesins, is about 
0.2 nm/ms if the supply of ATP is sufficient, and that 
v decreases with the lowering of ATP concentration fol- 
lowing a Michaelis-Menten type relation. 



FIG. 8: The stationary density profiles for loh = 0.1 (left) and 
ujh = 0.2 (right) in the case u a = 0.001. The blue and red 
lines correspond to the densities of state 1 and 2, respectively. 
The dashed lines are the mean-field predictions for periodic 
systems with the same parameters. 

In contrast to the phase dia gram s in the a — /3-plane re- 
ported by earlier investigators J48[ I4QL IHH ] . we have drawn 
the phase diagram of our model in the u a — UJh plane by 
carrying out extensive computer simulations for realis- 
tic parameter values of the model with open boundary 
conditions |59ll67|. The phase diagram shows the strong 
influence of hydrolysis on the spatial distribution of the 
motors along the MT. For very low u>h no kinesins can ex- 
ist in state 2; the kinesins, all of which are in state 1, are 
distributed rather homogeneously over the entire system. 
However, immobile (but fluctuating) shock is observed in 
the densiti profiles of kinesins if u>h is finite. 
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FIG. 9: Schematic description of the traffic of ribosomes 
where each ribosome has been treated as an "extended" par- 
ticle (£ = 2 in this figure. 



VI. RIBOSOME TRAFFIC AND PROTEIN 
SYNTHESIS 

Translation, the process of synthesis of proteins j6^ 
by decoding genetic information stored in the mRNA, 
is carried out by ribosomes. Understanding the physi- 
cal principles underlying the mechanism of operation of 
this complex macromolecular machine |69j | will not only 
provide insight into the regulation and control of protein 
synthesis, but may also find therapeutic a ppl ications as 
ribosome is the target of many antibiotics |7f| . 

Most often many ribosomes move simultaneously on 
the same mRNA strand while each synthesises a protein. 
In all the earlier models of collective traffic-like move- 
ments of ribosomes 

EH El El El El El E3 . 

the entire 

ribosome is modelled as a single "self-propelled particle" 
ignoring its molecular composition and architecture. The 
typical size of an individual ribosome is much larger than 
that of a single codon (i.e., a triplet of nucleotides). This 
feature of the ribosome is captured in most of the recent 
theoretical models by postulating that the size of the self- 
propelled particle is £ (£ > 1) where the unit of length is 
set by the size of a codon. Moreover, in these models the 
inter-ribosome interactions are captured through hard- 
core mutual exclusion and the dynamics of the system is 
formulated in terms of rules that are essentially straight- 
forward extensions of the TASEP [55j]. The qualitative 
features of the flow properties and spatio-temporal or- 
ganisation in these models are very similar to those of 
TASEP which corresponds to the special limit I = 1 . 

In reality, a ribosome has a complex architecture. Each 
ribosome consists of two subunits and has four binding 
sites. Of these, three sites (called E, P, A), which are 
located in the larger subunit, bind to tRNA, while the 
fourth binding site, which is located on the smaller sub- 
unit, binds to the mRNA strand. The translocation of 
the smaller subunit of each ribosome on the mRNA track 
is coupled to the biochemical processes (see fig llOj l occur- 
ing in the larger subunit. 

Let us begin the biochemical cycle with state 1 where 
a tRNA is bound to the site P of the ribosome. A tRNA- 
EF-Tu complex now binds to site A and the system makes 
transition from the state 1 to the state 2. As long as the 
EF-Tu is attached to the tRNA, codon-antiodon binding 
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Phosphate 



FIG. 10: A simplified biochemical cycle of s ribosome. 



can take place, but the peptide bond formation is pre- 
vented. The EF-Tu has a GTP part which is hydrolizcd 
to GDP, driving the transition from state 2 to 3. Fol- 
lowing this, a phosphate group leaves, resulting in the 
intermediate state 4. This hydrolysis, finally, releases 
the EF-Tu, and the peptide formation is now possible. 
After this step, EF-G, in the GTP bound form, comes 
in contact with the ribosome. This causes the tRNAs 
to shift from site P to E and from site A to P, site A 
being occupied by the EF-G, resulting in the state 5. 
Hydrolysis of the GTP to GDP then releases the EF-G 
followed by conformatinal changes that are responsible 
for pulling the mRNA-binding smaller subunit by one 
step forward. Finally, the tRNA on site A is released, 
resulting in completion of one biochemical cycle; in the 
process the ribosome moves forward by one codon 0. 

Very recently we EH have developed a quantitative 
model that not only incorporates the hard-core mutual 
exclusion of the interacting ribosomes, but also cap- 
tures explicitly the essential steps in the biochemical cy- 
cle of each ribosome, including GTP (guanine triphos- 
phate) hydrolysis, and couples it to its mechanical move- 
ment during protein synthesis. The modelling strategy 
adopted in ref . |78j for incorporating the biochemical cy- 
cle of ribosomes is similar to that followed in our earlier 
work [H^ on single-headed kinesin motors KIF1A. How- 
ever, the implementation of the strategy is more difficult 
in case of ribosome traffic not only because of the higher 
complexity of composition, architecture and mechano- 
chemical processes of the ribosomal machinery but also 
because of the heterogeneity of the mRNA track E3 ■ The 
details of our work on ribosome traffic has been reported 
elsewhere 17811. 
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VII. SUMMARY AND CONCLUSION 

In this paper we have presented examples of collective 
behaviours of molecular motors that emerge from the co- 
ordinations, cooperations and competitions at different 
levels in the sub-cellular world. In our original works 
so far we have synthesized the single-motor mechano- 
chemistry and multi-motor interactions to develope the- 
oretical models that make experimentally testable quan- 
titative predictions. In particular, we have developed 
models of KIF1A traffic on a MT track and ribosomc 



traffic on a mRNA track. 
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